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Abstract 

Although repulsive effects have been predicted for quantum vacuum forces between bodies with 
nontrivial electromagnetic properties, such as between a perfect electric conductor and a perfect 
magnetic conductor, realistic repulsion seems difficult to achieve. Repulsion is possible if the 
medium between the bodies has a permittivity in value intermediate to those of the two bodies, but 
this may not be a useful configuration. Here, inspired by recent numerical work, we initiate analytic 
calculations of the Casimir-Polder interaction between an atom with anisotropic polarizability and 
a plate with an aperture. In particular, for a semi-infinite plate, and, more generally, for a wedge, 
the problem is exactly solvable, and for sufficiently large anisotropy, Casimir-Polder repulsion is 
indeed possible, in agreement with the previous numerical studies. In order to achieve repulsion, 
what is needed is a sufficiently sharp edge (not so very sharp, in fact) so that the directions of 
polarizability of the conductor and the atom are roughly normal to each other. The machinery for 
carrying out the calculation with a finite aperture is presented. As a motivation for the quantum 
calculation, we carry out the corresponding classical analysis for the force between a dipole and a 
metallic sheet with a circular aperture, when the dipole is on the symmetry axis and oriented in 
the same direction. 
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I. INTRODUCTION 



There has been increasing interest in utihzing the quantum vacuum force or the Casimir 

n 

effect in nanotechnology employing mesoscopic objects Although the original Casimir 



effect, between parallel conducting or dielectric plates separated by vacuum [2|, |3|, always 
gives an attractive force between the plates, introducing a material (liquid) with an inter- 
mediate value of the dielectric constant can result in repulsion [4], which has now been 
observed |5.|]. For precursors, see {g-IO]. [The first experimental test of the Lifshitz theory 
with an intermediate liquid (helium) was that of Sabisky and Anderson [11]; application of 
the Lifshitz theory to the melting of water ice was considered by Elbaum and Schick 12|.] 
A recent experiment involving air bubbles in a liquid with boundary walls is described in 
Ref. [l3|. However, this type of repulsion is unlikely to have many applications in building 



devices. 

There are well-known repulsive quantum forces in vacuum. The first example was found 
by Boyer [l^. He computed the self-stress of a perfectly conducting spherical shell due to 
quantum electrodynamic field fluctuations and found a repulsive result, but the meaning of 
such a self-energy is extremely obscure. He later found 15(| a more observable effect, that the 
force between a perfect electrically conducting plane (e, the permittivity, goes to infinity) 
and a parallel perfect magnetic conducting plane {fi, the permeability, goes to infinity) is 
repulsive. This, again, may be a difficult situation to approximately replicate in practice, 
because the unusual magnetic properties must persist over a wide frequency range. 

There has been extensive interest in designing metamaterials that could give rise to 



Casimir repulsion by simulating a magnetic response |16l-l2l|. Despite some early optimism, 
the conclusion seems to have transpired that repulsion is impossible between metamaterials 
made from dielectric and metallic components 22|-|24| . For recent attempts using dielec- 
tric/magnetic setups see Refs. |25l-l27l|. 

Several years ago there was an interesting suggestion by Sopova and Ford 28| that the 
force between a small dielectric sphere and a dielectric wall was oscillatory, so there were a 
number of repulsive regimes. However, this effect was canceled by plasmon modes leaving 
the usual attractive result ^9|. Earlier Ford had suggested |30| that the frequency response 



of materials might be manipulated in order to achieve repulsion, but this was proved to be 
impossible 31 1. 
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Thus it was extremely interesting when Levin et al. showed examples of repulsion between 
conducting objects, in particular between an elongated cylinder above a conducting plane 
with a circular aperture 32|. (An analytic counterpart is given in [33'].) They first gave 
examples of repulsive forces between arrays of electric dipoles, and an electric dipole and a 
conducting plane with an aperture cut out. Then they turned to quantum vacuum forces 
between conducting objects, computed by quite impressive "brute force" finite-difference 
time-domain and boundary-element methods. 

The purpose of the present paper is to try to understand these phenomena analytically. 
We first show, in Sec. [TTl that there is no repulsion possible in the weak coupling regime, 
where because the materials are dilute one may sum Casimir-Polder interactions between 
atoms 34|. However, there is repulsion in classical electrostatics between a system of three 
dipoles (Sec. IIIII) and between a fixed dipole and a conducting plane with an aperture, which 
we discuss in Sec. HVt both in two and three dimensions. This is an interesting pedagogical 
problem, for it involves mixed coupled integral equations, like those for an electrified disk, or 
a plane with an aperture with different constant electric fields at large distances above and 
below the punctured plane |35|. These problems exhibit closed form solutions, and clearly 
exhibit repulsion when the dipole is directly above the aperture and is sufficiently close. In 
Sec. |V]we turn to the real problem, that of the Casimir-Polder force between an anisotropic 
polarizable atom and a punctured dielectric plane. Because solving the integral equations 
arising for the Green's dyadic for the plate with aperture is rather complicated, in Sec. I VI I 
we content ourselves with computing the Casimir-Polder interaction between a polarizable 
atom and a perfectly conducting wedge. When the opening angle of the wedge approaches 
2tt, this describes the interaction between an atom and a semi-infinite conducting plane. 
We exhibit situations in which repulsive forces in certain directions can arise for anisotropic 
atoms, in qualitative agreement with numerical work {32 1. In Appendix A we give another 
derivation of the Casimir-Polder energy formula for the wedge, based on a closed form for the 
Green's dyadic, and in Appendix B we give a classical calculation of a conducting ellipsoid 
above a conducting plate with a circular aperture in the presence of a background field. 

A word about terminology: When we say "atom" we mean any microscopic particle which 
may be described by a polarisability tensor. Our calculations assume that we are in the 
retarded regime, so that static (frequency-independent) polarizabilities may be employed. 
Should lower frequency transitions dominate (which could occur with some molecules), so 
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FIG. 1: Two-dimensional geometry of a needle of length L a distance Z above a line with a gap of 
width a. 

that the separations are in the non-retarded regime, electrostatic results are valid (but for 



a factor of 1/2— See Eq. f Hl5|) below and Ref. |36|). 

In this paper we set /i = c = 1, and all results are expressed in Gaussian units except 
that Heaviside-Lorentz units are used for Green's dyadics. 

II. WEAK COUPLING CALCULATION 
A. Scalar field 

We first illustrate the ideas by considering the case of a massless scalar field in two 
dimensions. The quantum vacuum energy between two weakly coupled potentials Vi and V2 
is 

the scalar analog of the Casimir- Polder force between atoms. Here we consider the potentials 
as shown in Fig. [H which represents a needle of length L on the symmetry axis a distance 
Z above a line with a gap of width a. The potentials are given by 

Vi{x,z) = Xi6{x)e{z - Z + L/2)e{Z + L/2 - z), (2.2a) 
V2{x,z) = X26{z)[e{x-a/2) + e{-x~a/2)]. (2.2b) 

This means that the interaction energy is 



rZ+L/2 




/ dz< 




1 Z-L/2 


./a/2 J -00 J 
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FIG. 2: Three-dimensional geometry of a polarizable atom a distance Z above a dielectric slab 
with a circular aperture of radius a. 



To get the force on the needle, we simply have to integrate on x, and differentiate with 
respect to the limits of the z integral: 



A1A2 



arctan(2Z/ a + L/a) arctan(2Z/ a — L/a) 



(2.4) 



2Z/a + L/a 2Z/a-L/a 

which, because F < 0, always represents an attractive force between the punctured line 
and the needle. Note that although the force vanishes at Z = 0, the energy there, which 
represents the work done in bringing the needle in from infinity, is not zero. 



B. Electromagnetic field 

Now we consider the quantum vacuum force between dilute dielectric media, which may 
be obtained from the Casimir- Polder potential between isotropic polarizable atoms , 



CP 



23 1 

Air r' 



(2.5) 



where r is the distance between the atoms. We might mention that equation f l2.5p is 
in general valid in the retarded limit where the atomic polarizability can be regarded as 
constant. (For more details, see the review 37|.) The result is applicable provided that 
the atom-plate separation is much greater than the atomic transition wavelength (typically 
some hundreds of nanometers for ground-state atoms). The media have dielectric constants 
Ei = l + AnNiai, where Ni represents the density of atoms of type i. Specifically, we consider 
a three-dimensional configuration, in which an atom of isotropic polarizability a is placed 
on the symmetry axis a distance Z above a dielectric plate of thickness t with a circular hole 
in the middle, as shown in Fig. [2l The quantum interaction energy is 
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(47r)2 
23 



tMS 



(dr) 



{t + 2Z){6a^ + {t + 2Zf) 



+ iZ ^ -Z) 



(2.6) 



eOvra^"'^ ^' [ (4a2 + (t + 2Z)2)3/2 
It is easy to see that the force F = —dU/dZ is always negative, i.e., attractive. 

A more favorable case for possible repulsion would be an anisotropic atom. It is easy to 
derive the appropriate generalization of the Casimir-Polder potential in this case, starting 
from the weak-coupling multiple scattering formula 



U 



12 



where the free Green's dyadic is {( = —iu) 



ro(r,r') = (VV-lC 



e-|CI|r-r'| 



(2.7) 



(2.8) 



Following the procedure given in Ref. 
anisotropic atom 
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47r|r — r'l 

we find for an isotropic medium facing an 



U 



1 



327r2 



(dr) 



slab 



R|7 



13trQ; + 7 



- R) ■ a ■ (r - R) 



(2.9) 



(r - R)2 

where R = (0, 0, Z) is the position of the atom, relative to the center of the aperture. This 
may be easily checked to reduce to the usual Casimir-Polder result (12. 5p when ex. = al. 

Let's consider the extreme case when only a^z is significant. Then the integrals may be 
easily carried out, with the result 



U 



eOvra^ 



t + 2Z 



[156a^ + 70a\t + 2Zf + 7(t + 2ZY] + {Z ^ -Z) 



(2.10) 



This, again, always gives rise to an attractive force. 

An interesting special case is when the aperture is small compared to the thickness of the 
dielectric. Then the energy is a step function, 
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U 



-azz{e-l)e{t-2\Z\), a^t, 



(2.11) 



307ra4 

which gives rise to a 5-function force just when the atom enters and exits the aperture. If 
the aperture is very large compared to the thickness of the slab, t ^ a, the energy and force 
are proportional to the thickness of the slab. 



U 



807ra4 



1 , ^ , 13a2 + 18^2 

a,,A£ - 1 



{a? + Z2)7/2 



aH. 



(2.12) 



7 




FIG. 3: Configuration of tliree dipoles, two of which are antiparallel, and one perpendicular to the 
other two. 

III. CLASSICAL DIPOLE INTERACTION 

It is possible to achieve a repulsive force between a configuration of fixed dipoles. Consider 
the situation illustrated in Fig. [31 Here we have two dipoles, of strength d2 and ^3 lying 
along the x axis, separated by a distance a. A third dipole of strength di lies along the z 
axis. If the two parallel dipoles are oppositely directed and of equal strength, 

da = -dg = (i2X, (3.1) 

and are equally distant from the z axis, and the dipole on the z axis is directed along that 
axis, 

di = diz, (3.2) 
the force on that dipole is along the z axis: 

F, - 3adid2 + ^2/4)7/2' (3.3) 

which changes sign at Z = a/4. That is, for distances Z larger than this, the force is 
attractive (in the —z direction) while for shorter distances the force is repulsive (in the +z 
direction). Evidently, by symmetry, the dipole-dipole energy vanishes at ^ = 0. Consistent 
with Earnshaw's theorem, the point where the force vanishes is an unstable point with 
respect to deviations in the x direction. 

In view of this self-evident finding, it might seem surprising that the interaction between 
a polarizable atom and a dilute medium (made up of polarizable atoms) studied in Sec. IIIBI 
failed to exhibit a repulsive regime, but this is because the medium is isotropic. 
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IV. CLASSICAL INTERACTION BETWEEN A DIPOLE AND A CONDUCTING 
PLANE WITH AN APERTURE 



In this section, we consider the interaction between a dipole and a perfectly conducting 
plane containing an aperture. We first consider two dimensions. (As above, we denote the 
Cartesian coordinates by x and z for uniformity with the three-dimensional situation.) 

A. Dipole above aperture in a conducting line 

If we use the Green's function which vanishes on the entire line ^ = 0, 



^(r, r') = - ln[(a; - x'Y + {z - z'Y] + ln[{x - x'Y + {z + z^ 



(4.1) 



so 



G{x,0;x',z') =0, (4.2) 
we can calculate the electrostatic potential at any point above the z = plane to be 



0(r)= / {dr')G{r,r')piT') + ^ [ dS'-^G{r,r'] 
Jz>o 477 Jap dz' 

where the volume integral is over the charge density of the dipole. 



(4.3) 



z'=0 



p(r) = -d- V(5(r-R), R=(0,Z). 



(4.4) 



The surface integral extends only over the aperture because the potential vanishes on the 
conducting sheet. If we choose d to point along the z axis we easily find 



{x,z> 0) = 2d 



z-Z 



z + Z 



x^ + {z-Zy x^ + {z + Z)\ 



1 



TC 



a/2 



dx'- 



a/2 



(x — x'Y + -2^ 



0(x',O), 



(4.5) 



where a is the width of the aperture. 

Now the free Green's function in two dimensions is 



Then the surface integral in Eq. (14.51) is 



oo I kx I 



(4.6) 



/TT 



(4.7) 
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in terms of the Fourier transform of the field 



dx e 



/ ^ CC 



'0(x',O) 



a/2 



dx' cos kxx'(f){x' , 0), 



(4i 



since 0) must be an even function for the geometry considered. Thus we conclude 

z — Z z + Z 



{x,z >0) = 2d 



_x-^ + {z-zy 



x^ 



[z + zY\ 



dk COS kxe-'''(j){k). (4.9) 



This becomes an identity as z — )■ 0. 
The electric field in the aperture is 



d 

E^{x, z = 0+) = -— z) 



Ad- 



X 



Z-0+ 



x^ + z^y 



vr 



dk k cos kx (f){k) . (4.10) 



On the other side of the aperture, there is no charge density, so for 2; < the potential is 



1 



x,z<0) = - / dk cos kxe (f){k), 



so the z-component of the electric field in the aperture is 



Ez{x, z = 0- 



d_ 

dz 



{x,z) 



1 



77 



dk k cos kx (j){k). 



(4.11) 



(4.12) 



2=0- "JO 

Because we require that the electric field be continuous in the aperture, and the potential 
vanish on the conductor, we obtain the two coupled integral equations for this problem. 



Ad- 



x^-Z^ _ 2 







dk k COS kx (j){k) , < |x| < a/2, 
dk cos kx (f){k) , \x\>a/2. 



In fact, these equations have a simple solution [39 1 

~ _ AZdTi Joikax/2) 

From this, we can work out the energy of the system from 



(4.13a) 
(4.13b) 

(4.14) 

(4.15) 



=Z,x=0 



where the factor of 1/2 comes from the fact that this must be the energy required to assemble 
the system. In computing this energy we must, of course, drop the self-energy of the dipole 
due to its own field. We are then left with 



d^ 



d 



int 



— / dkke-^^(t){k) 







4^2 27r 



4Z2 
'{o? + AZ^f 



1 







2 (x2 + 4Z2/a2)3 



(4.16) 
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where to get the second hne we used the derivative of Eq. f l4.17p . This is exactly two times 
larger that the result quoted in Ref. 32|.^ Since this vanishes at Z = and Z = oo, the 
force must change from attractive to repulsive, which happens at Z = a/2. 



B. Three dimensional aperture interacting with dipole 

It is quite straightforward to repeat the above calculation in three dimensions. Again 
we are considering a dipole, polarized on the symmetry axis, a distance Z above a circular 
aperture of radius a in a conducting plate. 

The free three-dimensional Green's function in cylindrical coordinates has the represen- 
tation 

r-, - = \ dkJo{kp)e~^\^\ (4.17) 

\J + Jo 

and so if we follow the above procedure we find for the potential above the plate 



0(r_L, z> Q) = d 



z — Z ^ z + Z 



\rl + {z- Zffl^ [r2 +(z + Z)2]3/2_ 

/■oo 

+ / dkke~^'Jo{kr^)^{k), (4.18) 
Jo 

where the Bessel transform of the potential in the aperture is 

/>oo 

$(fc)= / t/ppJo(M0(P,O). (4.19) 

Thus the integral equations resulting from the continuity of the z-component of the electric 
field in the aperture and the vanishing of the potential on the conductor are 

r2 - 2Z2 i-oo 



^ rJ2^^^2]5/2 = / dkk'^Jo{kr±)^k), r± < a, (4.20a) 







POO 

0= dkkJo{kr^)^{k), > a. (4.20b) 
Jo 



The solution to these equations is given in Titchmarsh's book '40'], and after a bit of 
manipulation we obtain 

- - (4.21) 



^ This is not the factor of 1/2 in Eq. ()4.15p . It is not possible to trace the origin of the discrepancy, since 
the authors of that reference merely quote the result. 
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Then the energy (14.150 may be easily evaluated using 



/ ..P/=e-.,,(*..) ^ (4.22) 



The energy can again be expressed in closed form: 



U 



arctan — + 



a Zl + ^I'i^Zjaf - {Z/a) 



Z a (l + ZVa2)3 



(4.23) 



This is always negative, but vanishes at infinity and at zero: 



4 

Z ^ : U ^ <f — . (4.24) 

This means that for some value of Z ~ a the force changes from attractive to repulsive. 
Numerically, we find that the force changes sign aX Z = 0.742358a. 

The reason why the energy vanishes when the dipole is centered in the aperture is clear: 
Then the electric field lines are perpendicular to the conducting sheet on the surface, and 
the sheet could be removed without changing the field configuration. 

Our goal is to analytically find the quantum (Casimir) analog of this classical respulsion. 

V. STRONG COUPLING FORCE BETWEEN AN ATOM AND A PUNC- 
TURED PLANE DIELECTRIC 

Now we turn to the real problem. Our starting point is the general expression for the 
vacuum energy 

f/=^Trlnrro-\ (5.1) 

where T is the full Green's dyadic for the problem, and Fq ^ is the inverse of the free Green's 
dyadic (12. 8p . namely 

To 1 = X V X -1. (5.2) 
In the presence of a potential V, the full Green's dyadic has the symbolic form 

r = (i-roV)-iro. (5.3) 

Here we are thinking of the interaction between a dielectric medium, characterized by an 
isotropic permittivity, so Vi = e — 1, and a polarizable atom, represented by a polarizability 
dyadic, as shown in Fig. [21 

V2 = 47ra(5(r-R), (5.4) 
12 



where R is the position of the dipole. We are only interested in a single interaction with 
the latter potential, so we have for the interaction energy 



U 



12 



TrV, 



■SVi 



--\n{l-ToV,) 



-Tr(ri-ro)V2, 



(5.5) 



where we have used Eq. (15.31) for the potential Vi describing the dielectric slab plus aperture 
and we have subtracted the term that represents the self-energy of the atom with its own 
field. This subtraction happens automatically if we start from the "TGTG" form, 

Uu = Trln (1 - riViFsVa) ^ Tr riViFoVa = '-Ti {T, - Tq) V2, (5.6) 

because V2 is weak. This implies the Casimir-Polder expression for the interaction between 
the polarizable atom and the dielectric 



CP 



/oo 
rfCtra-(r-ro)(Z, Z). 
'OO 



(5.7) 



(5.8) 



We could also derive this result from the formula for the force on a dielectric body in an 
inhomogeneous electric field |36|, 

Stt 

which classically says that a dielectric body experiences a force pushing it into the region of 
stronger field. This implies the interaction energy 

1 



U 



'-aE'{0,Z), 



and when we make the quantum- field-theoretic replacement. 



-^(E(r)E(rO)->ir(r,r') ^ 

An t I 



Ztt 



(5.9) 



(5.10) 



we recover the static isotropic version of Eq. (15. 7p after the self-energy is subtracted. 



A. No aperture 

When the aperture is not present, we are considering the well-studied case of a dielectric 
slab, of thickness t, interacting with a polarizable atom. Because the Green's dyadic in 
this situation, denoted r'^°\ then possesses translational invariance in the x-y plane, we can 
express it in terms of a reduced Green's dyadic, 

r(°nr,r') = J ^^e^^-<^-^'^-g{z,z';k^). (5.11) 
13 



In the case of an isotropic atom, the trace of the Green's dyadic occurs, which is for the 
reduced Green's dyadic 



tr g(Z, Z) = -eg'^iZ, ^) + + ^i) 9^iZ, Z') 



(5.12) 



Z'=Z 



in terms of the transverse electric (H) and transverse magnetic (E) Green's functions. These 
subtracted quantities are for z, z' above the dielectric 

H,Ef^ i\ H,E/ '^ 1 



9 ' [z,z )-gf^' [z,z 
in terms of the reflection coefficients 



2k 



_j^H,E^-K{z+z'-t) 



H 



R 



K — K 



KK 



' 1 



K — r' , kR' 1 
— + 4t ^, 



(5.13) 

+ 4 (5.14a) 

(5.14b) 

(5.15) 

D^[^^,) e--l, (5.16) 

with D obtained from this by replacing n' by R' except in the exponent. These results are 
rather trivially obtained by multiple scattering arguments. 
Now the interaction energy is 



where 



and 



K + k' 



U = -a 



dC 



a 
An 



oo 
oo 



dC, I dki—e 



+ '-^(2e + e)[k-'[i + -]+c 



(k + k' /eYe^'^'^ — {k — k' /ey 



This is precisely the result found, for example, by Zhou and Spruch 



.(5.17) 



B. Integral equations for Green's dyadic 

We now specialize to the case where the plane z = consists of a perfectly conducting 
screen with a circular aperture of radius a at the origin. The Green's dyadic satisfies the 
differential equation 

X V X -1 ) • r(r - r ) = l(5(r - r'), (5.18) 



a;2 
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subject to the boundary conditions 

z X rrr,r'^ 



0, 



(5.19) 



rx|>a,2=0 

which just states that the tangential components of the electric field must vanish on the 



conductor. Following Levine and Schwinger 42| we introduce auxiliary electric and magnetic 
Green's dyadics r^^'^''(r, r') which satisfy the same differential equation (15.181) but with the 
boundary conditions satisfied on the entire z = plane: 



z X r(i)(r,r') 



0, z X (V X r(2)(r,r')) 



z=0 



0. 



(5.20) 



2=0 



These can be constructed in terms of the free Green's dyadic Fq, subject only to outgoing 
boundary conditions at infinity, as given in Eq. (12.81) . 



ro(r,r') = (lu;2 + VV) G(|r-r'|), 
expressed in turn in terms of the Helmholtz Green's function 



GiR) 



^i\ui\R 



(5.21) 



(5.22) 



We can write, after the Euclidean rotation \uj\ — )■ iC,, the free Green's dyadic in the explicit 
form (R = r — r') 



rn(r,r') 



G{R) 



l{l + CR + eR') 



RR 



(3 + 3Ci? + C^^^) 



(5.23) 



In terms of this last dyadic, the auxilliary Green's dyadics have the form 

z,z' >Q: r«'(2)(r,r') = r(°)(r, r') T T^^Hr, r' - 2z^') • (1 - 2zz). (5.24) 



Now using Green's second identity, it is easy to prove 

V X r(2)(r,r') = [V X r«]^(r',r). 



r«-(2)( 



r, r 



[r(i).(2)irfr',r) 



(5.25a) 
(5.25b) 



where T signifies transposition. In the same way we may derive the integral equations for 
the Green's dyadic for the screen with the aperture 



^, > : r(r, r') = T^^) (r, r') - ^ j dS"V x T (r, r") • z x r(r", r 



z<0<z': r(r,r') 



^ I dS"V X rL'^(r, r") • z x r(r", r'). 



(5.26a) 
(5.26b) 
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where the ± subscripts on T^'^^ indicate that the Green's function is defined in the domain 
above or below the z = plane. The continuity of the z-component of the electric field in 
the aperture then leads to the integral equation 



r(i), 



r, r 



^ I dS"z ■ V X (rj^ + rL'^)(r, r") ■ z X r(r", r'). (5.27) 



z^0+ "5 ^ ap 

The system of integral equations defining the Green's dyadic is rather more complicated 
that those describing the corresponding (classical) static potential problem, so we will defer 
the discussion of strategies for its solution to a subsequent publication. We will here turn 
to a situation that can be solved exactly. 

VI. CASIMIR-POLDER FORCE BETWEEN ATOM AND A CONDUCTING 
WEDGE 

The interaction between a polarizable atom and a perfectly conducting half-plane is a 
special case of the vacuum interaction between such an atom and a conducting wedge. For 



the case of an isotropic atom, this was considered by Brevik, Lygren, and Marachevsky 



43|. 



(This followed on earlier work by Brevik and Lygren 4^ and DeRaad and Milton 45|.) 
In terms of the opening dihedral angle of the wedge f2, which we describe in terms of the 
variable p = tt/Q, the electromagnetic Green's dyadic has the form (here the translational 
direction is denoted by y, and one plane of the wedge lies in the z = plane, the other 
intersecting the xz plane on the line 6 = Q — see Fig. Hj) 



r(r.r') = 2pV' / f [ - .#.#"(Vi - t^)^F„,(p.p-) '=°^'"P'"'°""P'' e'>l-') 

m=0 L 

(6.1) 



1^ „/Nsinmp^sinmp^'^.fc(^_^,)' 



+ ^^'*-Gmp{p,p') —e 

00 TT 

The first term here refers to TE (H) modes, the second to TM (E) modes. The prime on 
the summation sign means that the m = term is counted with half weight. In the polar 
coordinates in the xz plane, p and 9, the H and E mode operators are 

^ = (6.2a) 
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FIG. 4: Polarizable atom, located at polar coordinates p, 0, within a conducting wedge with 
dihedral angle $7. 

where the transverse Laplacian is 

In this situation, the boundaries are entirely in planes of constant 6, so the radial Green's 
functions are equal to the free Green's function 

^F^,(p,p') = -Gmp{p,p') = -^Jmp{Xp<)H^l{Xp^), (6.4) 

CO CO A A 

with = co'^ — k'^ . We will immediately make the Euclidean rotation, co ^ iC, where A — >■ in, 
k"^ = C"^ + A;^, so the free Green's functions become —K,^^Imp{Kp<)Kmp{i^p>)- 

We start by considering the most favorable case for CP repulsion, where the atom is 
only polarizable in the z direction, that is, only a^z 7^ 0. In the static limit, then the only 
component of the Green's dyadic that contributes is 

Ti-^zz = -r^ / dkd(< \('^ sm^ 9 sm^ mp9 - k'^ cos^ 9 cos^ mp9] Impi'^P<)Kmp(np>) 

27: An-^ J l"- k^P<P> 

— [k^ sin^ 9 sin^ mp9 — cos^ 9 cos^ mp9] I'j^p{Kp<:)K!^p{Kpy)^ . (6.5) 

Here we note that the off diagonal p-9 terms in T cancel. We have regulated the result by 
point-splitting in the radial coordinate. At the end of the calculation, the limit p< ^ p> = p 
is to be taken. 

Now the integral over the Bessel functions is given by 



/ dKKh{np<)K^{hip>) ^ — — , (6.6) 
Jo - ? ) 
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where ^ = p^jp-^. After that the m sum is easily carried out by summing a geometrical 
series. Care must also be taken with the m = term in the cosine series. The result of a 
straightforward calculation leads to 



cos Id 



+ finite, 



(6.7) 



where the divergent term, as ^ — )■ 1, may, through a similar calculation, be shown to be that 
corresponding to the vacuum in absence of the wedge, that is, that obtained from the free 
Green's dyadic. Therefore, we must subtract this term off, leaving for the static Casimir 
energy (15. 7p 



JJZZ 



«..(0) 



Stt p4 sin p6' 



V 



— 1) sin^p6' + 



(p^-l)(p^ + ll) 
45 



sin"^ "pQ cos 26' 



(6.8) 



This result is derived by another method in Appendix lAl 

A small check of this result is that as 6* — )■ (or — )■ 1]) we recover the expected Casimir- 
Polder result for an atom above an infinite plane: 

a..(0) 



TTZZ , 

Uqp — ^ 



(6.9) 



in terms of the distance of the atom above the plane, Z = p9. This limit is also obtained 
when p — )■ 1, for when = tt we are describing a perfectly conducting infinite plane. 

A very similar calculation gives the result for an isotropic atom, a = al, which was first 



given in Ref. 

[/cP = 



43|: 



3a(0) 



P 



(6.10) 



Svrp^ sin^ pO 

Note that this is not three times f/^p in Eq. (16. 8 p because the cos 26' factor in the last term 
in the latter is replaced by —1/3 here. This case was reconsidered recently, for example, in 
Ref. y. 



A. Repulsion by a conducting half-plane 



Let us consider the special case p = 1/2, that is = 27r, the case of a semi- infinite 
conducting plane. This was the situation considered, for anisotropic atoms, in recent papers 
by Eberlein and Zietal 47H49| . Note that in such a case, for the completely anisotropic 
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FIG. 5: Polarizable atom, above a half conducting plane, free to move on a line perpendicular to 
the plane but a distance X to the left of the plane. 

atom, [/qp = at ^ = 7r/2, that is, there is no force on the dipole when it is polarized 



48 



49|. 



perpendicular to the half-sheet and directly above the edge, as observed in Refs. 

Consider a particle free to move along a line parallel to the distance X to the 

left of the semi-infinite plane. See Fig. [51 The half-plane x < constitutes an aperture of 
infinite width. With X fixed, we can describe the trajectory by m = X/ p = — cos 6, which 
variable ranges from zero to one. The polar angle is given by 

e l + u 



sin^ - 



2 2 

The energy for an isotropic atom is given by 

a(0) 1 



(6.11) 



CP 



327r X* 



V{u) 



where 



V{u) = 3u^ 



1 



[l + u) 



u 



The energy for the completely anisotropic atom is 

1 

If we consider instead a cylindrically symmetric polarizable atom in which 



(6.12) 



(6.13) 



(6.14) 



a. = a^^zz + 'ja^^{xk + yy) = a^^{l - 7)zz + 'ja^^l, 



(6.15) 



where 7 is the ratio of the transverse polarizability to the longitudinal polarizability of the 
atom. Then the effective potential is 



(i-7)v;. + 7^, 
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(6.16) 



and the z-component of the force on the atom is 



327r X5 



d_ 

du 



V(l-7)(l-3«^) + i(l + 27)n«) 



(6.17) 



where V is given by Eq. f l6.13p . Note that the energy f l6.16p . or the quantity in square 
brackets in Eq. fl6.17p . only vanishes at m = 1 (the plane of the conductor) when 7 = 0. 
Thus, the argument given in Ref. 32] applies only for the completely anisotropic case. 

The force is plotted in Figs. El [71 It will be seen that if 7 is sufficiently small, when 
the atom is sufficiently close to the plane of the plate the z-component of the force is 
repulsive rather than attractive. The critical value of 7 is 7c = 1/4. This is a completely 
analytic exact analog of the numerical calculations shown in Ref. |3^], where the interaction 
was considered between a conducting plane with an aperture (circular hole or slit), and a 
conducting cylindrical or ellipsoidal object. Our calculation demonstrates that three-body 
effects are not required to exhibit Casimir-Polder repulsion. 

It is interesting to observe that the same critical value of 7 occurs for the nonretarded 
regime of a circular aperture, as follows from a simple computation based on the result of 



Ref. [495. For example, applying the result there for an atom with polarizability given by 
Eq. (16.151) placed a distance Z along the symmetry axis of an circular aperture of radius a 
in a conducting plane gives an energy 



U 



167r2 



dC, a^^(C) 



X 



Z3 



(I + 7) 



TT 



arctan ■ 



2aZ 



2aZ 



{Z^ + a?f 



;i+7)(^^-a^ 



'l-l)a'Z' 



(6.18) 



It is easy to see that this has a minimum for 2; > 0, and hence there is a repulsive force close 
to the aperture, provided 7 < 7c = 1/4. 



B. Repulsion by a wedge 

It is very easy to generalize the above result for a wedge, p > 1/2. That is, we want 
to consider a strongly anisotropic atom, with only a^z significant, to the left of a wedge of 
opening angle 

/3 = 27r-fi, (6.19) 
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FIG. 6: (Color online) The z-component of the force between an anisotropic atom (with ratio of 
transverse to longitudinal polarizabilities 7) and a semi-infinite perfectly conducting plane, z = 0, 
X > 0. Fz = — Q22/(327rX^)/(ti) in terms of the variable u = X/p = —cos9. Here the atom lies 
on the line x = 0, y = —X, and p is the distance from the edge of the plane and the atom. Here, 
/ > corresponds to an attractive force on the z direction, and / < corresponds to a repulsive 
force. The different curves correspond to different values of 7, 7 = to 1 by steps of 0.1, from 
bottom to top. For 7 < 1/4 a repulsive regime always occurs when the atom is sufficiently close to 
the plane of the conductor. 

as shown in Fig. [HI We want the z axis to be perpendicular to the symmetry axis of the 
wedge so the relation between the polar angle of the atom and the angle to the symmetry 
line is 



where, as before, 6 is the angle relative to the top surface of the wedge. Then, it is obvious 
that the formula for the Casimir-Polder energy (16. 8p is changed only by the replacement 
of cos 2^ by cos 20, with no change in sinp6'. Now we can ask how the region of repulsion 
depends on the wedge angle /3. 

Write for an atom on the line x = —X 



^ = 6 + /3/2 



(6.20) 



CP ~ 



'22 



^22(0) 

8ttX^ 



(6.21) 



21 




0.990 0.992 0.994 0.996 0.998 1.000 



u 

FIG. 7: (Color online) Same as Fig. [6l The region close to the plane, 1 > u > 0.99, with 7 near 
the critical value of 1/4. Here from bottom to top are shown the results for values of 7 from 0.245 
to 0.255 by steps of 0.001. 




FIG. 8: A polarizable atom outside a perfectly conducting wedge of interior angle /?. The atom is 
located at polar angles p, (p relative to the symmetry plane of the wedge. 



where 



V{4>) = cos^ 



p 



2 p^p^ - 1) 1 . . 

^ +— (/-I)(j9^ + ll)cos20 



.^"^ 2 TT-13/2 2 7r-/3/2 



(6.22) 



At the point of closest approach. 



V(7r) = -(4p2-l)(4p2 + ii), 



(6.23) 



so the potential vanishes at that point only for the half-plane case, p = 1/2. The force in 
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2.0 2.5 3.0 

FIG. 9: (Color online) The z-component of the force on an completely anisotropic atom moving 
on a line perpendicular to a wedge. The different curves are for various values of (3 from to vr by 
steps of 7r/20, from bottom up. The last few values of (3 have a markedly different character from 
the others. 



the z direction is 

/(</.) = cosV^^. (6.24b) 

Fig. |9] shows the force as a function of for fixed X. It will be seen that the force has a 
repulsive region for angles close enough to the apex of the wedge, provided that the wedge 
angle is not too large. The critical wedge angle is actually rather large, /3c = 1.87795, or 
about 108°. For larger angles, the 2;-component of the force exhibits only attraction. Of 
course, the force is zero for /3 = vr because then the geometry is translationally invariant in 
the z direction. 



VII. CONCLUSIONS 



This paper may be thought of as a counterpart to Ref. 32|. While that reference pro- 
ceeded on the basis of numerical calculations, we have used analytic approaches. After 
some examples indicating that Casimir-Polder attraction is typical, and always seems to 
occur in weak coupling, we demonstrate that the quantum-vacuum Casimir-Polder interac- 
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tion for a sufficiently anisotropic atom above a conducting half plane can exhibit regimes of 
repulsive forces for motion confined to certain specified directions. This directly translates 
into repulsion between such an atom and a plane with an aperture for motion along a line 
perpendicular to the plane. More complete analysis of that case will be presented elsewhere. 

As we were putting finishing touches on this paper, Ref. 49| appeared, which demon- 
strates in the nonretarded (van der Waals) regime, repulsion could occur between an 
anisotropically polarizable atom and a conducting plate with an aperture. The critical 
value of the anisotropy is the same as found here. 

Perhaps most remarkable here is that not only can we achieve repulsion with a half-plane, 
but also with a wedge geometry, even when the interior angle of the wedge is greater than 90°. 
This indicates that while anisotropy in both the atom and the conductor must be present 
for repulsion, the anisotropy in the latter need not be too extreme, and that repulsion in 
other geometries may be readily achievable. Three-body forces are not required, nor is a 
high degree of symmetry, as was present in Refs. 32|, |49 |. 
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Appendix A: Derivation of anisotropic wedge CP force from closed-form Green's 
function 



Many years ago Lukosz gave a closed form for the Green's functions for a perfectly 
conducting wedge SOj. The four- dimensional Euclidean Green's dyadic has the closed form 



r(r -t\x- x',p,p',> 



MM'G" + NN'G 



(Al) 



where the transverse differential operators are [cf. Eq. (16. 2p ] 



M 



p d(j) 



op 



N- p— + ^1-^ 

dp p 90 ' 



(A2) 
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where there is an additional contribution to N in the x direction. This Green's dyadic is 
the frequency Fourier transform of that discussed in Sec. IVII Here the E (TM) and H (TE) 
Green's functions have the form 

for a wedge of dihedral angle Q, with G [—Q/2, Q/2]. Here 

y(x o t- x' d t'- tb) - ^ sinh(7rt;/l^) 

Q-Kilpp' smht; cosu.[Tiv — cos^Tr-i/^/iZ) 



where 



sinh — = - 
2 2 



{T-T'f + {X-X'f + {P-P'f 



1/2 



For the interaction with an atom possessing only an polarizability, we need 



(A5) 



+ 2 ( sin cos <P'\§^^ + sin 0' cos 0^^^ ) x(0 - 0') 



cos 



f I d d d d \ 



Here, we have suppressed all the arguments in x except for the angular ones. For our 
application here, we are interested in the coincidence limit, so from the outset we can set 
T = t' and X = x'. Then 

.nKl^l^. e^^, (AT, 

which implies 

i; = -ln^. (A8) 

Now we expand first in — 0', then after the differentiations set = 0', and then expand in 
V, that is, in 1 — We immediately note that the mixed derivative term in Eq. (IA6P does 
not contribute, because there is no linear term in — 0'. The result of a straightforward 
calculation is 

cos 2^ J 16 1 2 2 , TIN 



1 r 2p2(p2_l^ 

IGtt^P^ \sin'^]96' 3 sin^ 
where p = n/Q, and we have switched to the angle from the "upper" plate, 9 = (j) + Q/2, 
which is chosen to run from to Q. The first term in Eq. (IA9p corresponds to the x(0 — 0') 
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contribution, and the second to the cos(0 + (p' — Q) contribution. Note, the divergent term 
(as ^ — > 1) is precisely the vacuum term given in Eq. (16.71) . and should be subtracted off, 
and the rest, when multiplied by —27iazz, coincides with Eq. 



Appendix B: Electrostatic aspects: Conducting ellipsoid outside a conducting plate 
with a circular hole 



Consider a conducting uncharged solid ellipsoid with semiaxes c > a > b, centered at 
X = Y = Z = 0. The ellipsoid is orientated such that the major semiaxis c lies along the 
Z axis. To describe the electrostatic potential in the external region, one can make use of 
ellipsoidal coordinates ^, r], (, corresponding to solutions for u of the cubic equation 



^2 ^2 y2 



+ u + u + u 



1. 



(Bl) 



The coordinate intervals are 



(B2) 



The relationships between the ellipsoidal and the Cartesian coordinates are given in Ref. 51 1 
and will not be reproduced here. We shall however need the line element. 



where 



h 



2R 



2Rn 



2R 



Ir^/iC-m-v), Rl = iu + c')iu + a')iu + b'), 



(B3) 

(B4) 
(B5) 



with u = ^,r],C- 

In the following we assume axial symmetry around the Z axis. Then a ^ b, t] —6^, 
and the equation for the surface of the ellipsoid becomes 

Z^ i?2 



with i?^ = + We now have 



Z = ± 



(e + c^)(C + c^) ' 



+ 



1/2 
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(B6) 



R 



{^ + b'){C + b') 

b'^-c^ 



1/2 



(B7) 
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The ellipsoidal coordinates ^, 77, ( reduce in the case of axisymmetry to so-called prolate 
spheroidal coordinates ^ and (, lying in the intervals 

00 > ^ > -b\ -b^>C> -c\ (B8) 

Surfaces of constant C, and ( are prolate spheroids and hyperboloids of revolution, the surfaces 
intersecting orthogonally. On the Z axis {R = 0) one has ( = —6^, Z = ±a/^ + c^, whereas 



in the XY plane {Z = 0) one has ( = —c^,R = a/^ + 6^. On the surface of the ellipsoid, 
In free space outside the ellipsoid the Laplace equation reads 



R, d 



^ R,^ 



Rt d 



^ Rr^ 



0. 



(B9) 



Assume now that the ellipsoid is placed in an external potential (pQ, axisymmetric with 
respect to the Z axis so that 0o = 0o(^, 0- We write the resulting potential in the form 

0(e,c) = 0o(e,c)[i + i^(O], (BIO) 

so that (PqF is the perturbation of the external field. As the boundary condition ^ = on 
the surface has to hold for all values of C, it is natural to make the ansatz that F depends 
on ^ only. 

Inserting Eq. f lBlOp into Eq. f lBQp we find that the terms containing F as a factor sum 
up to zero, the reason being the validity of Eq. (IB9I) also when is replaced by (pQ. The 
remaining terms containing F'{C,) and F"{C,) yield the equation 



d^F dF d ^ n 



(Bll) 



When integrating this equation, in order to preserve the validity of the ansatz F = F{C,), 
the coordinate ( in 0o has to be regarded as a parameter. The integration thus has to 
extend from ^ = (the surface) in the outward direction, along a line on the hyperboloid C 
= constant. 

The solution of Eq. fIBlip can be written as 

dC 



A 



R^<pr 



(B12) 



where the constant A is determined from the condition -F(O) = —1 on the ellipsoid surface. 
That means. 



Jo 



(B13) 
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We now specify the form of 0o, as the potential from a grounded conducting plate lying in 
the xy plane, when far from the plate there are constant electric fields, directed normal to 
the plate, having different values on either side. In the plate there is a circular opening with 
radius a (this radius not to be confused with the semiaxis a mentioned above). The center 
of the opening is at position x = ?/ = z = 0. It is known (Ref. [35|, Sec. 3.13) that on the z 
axis 



1 



\z\ 
a 



arctan ■ 



(B14) 



(t)o{z) = $00 

where $00 is a constant. At the origin, 0o = $00 • At infinity, \z\ ^ 00, (po ^ 0. 

The center of the vertically oriented ellipsoid is at position z = zq. Thus z = zo + Z. We 
will assume that the ellipsoid is so slender that the variation of 0o in the transverse x and y 
directions can be neglected. Thus we adopt the expression ( ]B14I) in the external field region 
of interest, 0o = 0o(^) 0^ ^ and ( being restricted to the same intervals (IBSP as before. 

We consider now the upper half of the ellipsoid, z > zq ot Z > 0. The nonperturbed 
potential, called 0o+) is then 



00 



arctan ■ 



(B15) 



Thus the potential 0+ in Eq. (IB 130 can be found numerically, inserting 0o+ together with 
= (^ + b'^)\/^ + c^. [In practice the following expansion can here be useful 52 1 



1 

— arctan x 

X 



< X < 1, 



(B16) 



k=l 



with coefficients a2k of order unity or less.] 

The induced surface charge density 0"+ on the ellipsoid is 



'eo d(f)+' 




2eo&c 90+ 




[hi \ 


5=0 







(B17) 



since on the surface hi = {b/2R^)^/^ = {l/2bc)^/^. In view of the relationships between 
the ellipsoidal and Cartesian coordinates this can be reexpressed as 



-2eo 



'Z^ R^' 


-1/2 










[ \ 


C=o 



(B18) 



From Eq. (IB 13 1) it follows that the derivative [d(f)o+/dC,]^=o does not contribute to a+ (recall 
that F{0) = —1). The remaining term is 
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1 -1 



(B19) 



Thus for z > ZQwe get as solution 



c4 



-1/2 



(B20) 



where (To+ is the constant 



2eo ^ 
62 $00 



rex 
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(B21) 



(recall again that a is the radius of the hole). The dependence of on the coordinates 
Z and R in Eq. ( 1B20I) is actually the same as for a charged ellipsoid in free space 51 1. 
The surface force density on the ellipsoid is ((j2/2eo)n, n being the outward normal. The 
slope of the tangent to the surface is dZ/dR = —{c^/lP')R/Z; the slope of n is accordingly 
{\? j (?^Z I R. Denoting this as tan^, we get when going over to ellipsoidal coordinates, 



tan 6' 



C + c2 



1/2 



(B22) 



The component of n along the Z axis is then 



nz = sin 9 



X + c' 

-c 



1/2 



(B23) 



and we can now find the total vertical force Fz+ on the upper half of the ellipsoid by 
integrating over the actual surface. The line element along the meridian is 



hdC = ^ 



C 



1/2 



and the surface element dA becomes 

dA = 2ttR hi dC 



nb 



nl/2 



dC 



As cr+ in Eq. (IB20p can be reexpressed as 



(B24) 
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(B26) 



we can calculate Fz+ as 



-nz dA 



2en c2 



irb'^ c 



(B27) 



The expression is positive as expected; the force is acting upwards. The only dependence 
on the position zq lies in (To+, as ao+ = ao^{zo) according to Eq. (1B2ip . 
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0o(p,O) 



The lower half of the ellipsoid, Z < 0, can be treated in an analogous way. A complicating 
element is here the presence of the conducting plate in the xy plane, for radii p > a. It 
means that we can no longer extend the integration over ^ in the solution (lB12p to infinity 
in a straightforward way. We observe that the undisturbed potential in the xy plane can be 
written as 

' (B28) 

0, p > a, 

where = + y^, $00 being the potential at the center. 

Our approach will be based on the following two assumptions: 

1) The ^ integration will be terminated on the xy plane, this implying that the effect of 
the perturbation is assumed to be small at that level. This approximation is expected to be 
good except when the distance between the lower end of the ellipsoid and the plane is small. 

2) Secondly, the integration over ^ will be assumed to run over trajectories lying close 
to the z axis, corresponding to ( = —6^. This assumption simplifies the mathematical 
analysis. It is supported by physical considerations also, since when the ellipsoid is slender 
the hyperboloids ( = constant emerging from the surface of the ellipsoid near its lower end 
become concentrated in the vicinity of the z axis. 

As according to Eq. (1B7I) the plane position z = in general corresponds to 
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(B29) 



our approximations imply that the ^ integration is terminated at 



^pla 



(B30) 



i.e., the same constant for the whole lower half of the ellipsoid. 
As solution for the perturbed potential we thus get 



•Cpla 



where 
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The force Fz~ on the lower half can now be calculated. As before, i?^ = + b'^)^/^ + c^. 
Equation (IB26P becomes replaced by 

b 



a_ = (To-- 



where now 



CTq- 
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The total force on the ellipsoid becomes 

Fz = Fz+ + Fz- 
which can be rewritten as 
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In the limiting case of a sphere, b ^ c, the expression becomes somewhat simpler, 
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We have made some numerical checks of these expressions (using Maple). They indicate 
that there is no change in the sign of the force for various input parameters for the geom- 
etry. The force is attractive, as expected. It turns out that the dependence on the upper 



integration limit ^pia 



is weak, as anticipated above. 
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